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Abstract 

The complete understanding of the basic constituents of hadrons and the hadronic 
dynamics at high energies are two of the main chahenges for the theory of strong 
interactions. In particular, the existence of intrinsic heavy quark components in the 
hadron wave function must be confirmed (or disproved). In this paper we propose 
a new mechanism for the production of D-mesons at forward rapidities based on 
the Color Glass Condensate (CGC) formalism and demonstrate that the resulting 
transverse momentum spectra are strongly dependent on the behavior of the charm 
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distribution at large Bjorken x. Our results show clearly that the hypothesis of 
intrinsic charm can be tested in pp and p{d)A collisions at RHIC and LHC. 

Key words: Quantum Chromodynamics, Intrinsic Charm, Color Glass Condensate 
PACS: 12.38.-t, 13.38.Bx, 13.60.Le 



1 Introduction 

A complete knowledge of the partonic structure of the hadrons is fundamen- 
tal to make predictions for both the Standard Model and beyond Standard 
Model processes at hadron colliders. Since the early days of the parton model 
and of the first deep inelastic scattering (DIS) experiments, determining the 
precise form of the quark and gluon distributions of the nucleon has been a 
major goal of high energy hadron physics. Over the last 40 years enormous 
progress has been achieved. In particular, data from HERA have dramatically 
improved our knowledge about the small-x behavior of the parton distribu- 
tions functions (PDF's), finding a strong rise of the structure functions in 
this kinematical range, which implies a high partonic density in the proton. 
From theoretical considerations it is expected that in the high energy limit the 
small- a; gluons in a hadron wavefunction should form a Color Glass Conden- 
sate (CGC). The current experimental results at HERA and RHIC provide 
some evidence for the CGC physics (For a review see, e.g., [Ij). Although they 
are not yet conclusive, these experiments show the manifestation of non-linear 
(saturation) effects in high energy processes. The measured multiplicities and 
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rapidity distributions in heavy ion collisions, the high pt suppression of parti- 
cle production at large rapidities in deuteron-gold collisions, geometric scaling 
and other observations [1] point towards saturation. Of course, so far one has 
been able to explain these measurements with other approaches. For example, 
in Ref. [2j it was shown that geometric scaling can be derived from the DGLAP 
equations [3] . From our point of view, the most important consequence of these 
observations is that a change of formalism seems to be required. In particular 
it seems to be necessary to replace the collinear factorization formula by a 
more powerful formalism, able to accommodate saturation effects if and when 
they are present. 

Another important improvement has occurred in our knowledge about the 
heavy quark contribution to the proton structure. In the last years several 
groups have proposed different schemes to determine these distributions con- 
sidering that the heavy quark component in PDF's can be perturbatively 
generated by gluon splitting (See e.g. [1]). This component is usually denoted 
extrinsic heavy quark component. Moreover, the possibility of an intrinsic com- 
ponent has been studied in detail and included in the recent versions of the 
CTEQ parameterization [5]. The hypothesis of intrinsic heavy quarks (IHQ) 
is a natural consequence of the quantum fluctuations inherent to Quantum 
Chromodynamics (QCD) and amounts to assuming the existence of a QQ 
{Q = c, b, t) as a nonperturbative component in the hadron wave function. 
One of the most striking properties of an IHQ state, such as \uudQQ), is that 
the heavy constituents tend to carry the largest fraction of the momentum 
of the hadron. Consequently, in contrast to heavy quarks produced through 
usual perturbative QCD, which emerge with small longitudinal momentum, 
the intrinsic charm component gives rise to heavy mesons with large fractional 
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momenta relative to the beam particles. Therefore, the existence of an intrinsic 
component modifies, for instance, the xp and rapidity distribution of charmed 
particles (See, e.g. [B]). Moreover, it can also lead to Higgs production at high 
[7] . Many of these features have been discussed in the pioneering works on 
intrinsic charm p[T0|9] . 

Although the existence of IHQ fluctuations in the proton has substantial and 
growing experimental and theoretical support, as revised in Ref. [11] and 
briefly discussed below, more definite conclusions are still not possible. One 
of the shortcomings of previous studies is the use of the QCD factorization 
at large Feynman xp, where it is expected to breakdown (See, e.g., [12j). Ba- 
sically, at large Xp the kinematics is very asymmetric, with the hadrons in 
the final state emerging from collisions of projectile partons with large light 
cone momentum fractions {Xp — )■ 1) with target partons carrying a very small 
momentum fraction {xt ^ 1). In this case, we have a scattering of a dilute pro- 
jectile on a dense target, where the small- a; effects coming from CGC physics 
are expected to show up and the usual factorization formalism is expected to 
breakdown [13j. This kinematical region has been explored in dAu collisions 
at forward rapidities at RHIC [H] and the suppression of high px hadron 
yields, anticipated on the basis of CGC ideas, was observed. The satisfactory 
description of these experimental data [THTB] gives us a strong indication that 
the appropriate framework to calculate particle production at large xp (large 
rapidities) is one based on CGC physics. Having this in mind, in this paper 
we adapt the formalism proposed in [T7j (to calculate the production of high- 
X quarks in dAu collisions) to D meson production in pp {p{d)A) collisions 
at forward rapidities, including intrinsic charm quarks in the projectile wave 
function. In this approach the target is treated as a dense system while the 
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projectile proton is taken to be a dilute system in the spirit of standard QCD. 
We shall assume that the charm quarks are already in the projectile and their 
densities are those given by the CTEQ parameterization [5]. Moreover, we will 
describe the interaction of these charm quarks in terms of the dipole scatter- 
ing amplitude recently proposed in [15] to describe the RHIC data. Finally, 
we will assume that the charm quarks fragment into D mesons according to 
the fragmentation functions given in Ref. [18]. Our goal is to estimate the 
Z^-meson spectrum at forward rapidities at RHIC and LHC energies assuming 
the presence (or absence) of an intrinsic charm component and evaluate how 
much this quantity is modified. We emphasize that in our calculation there is 
no new free parameter. 



In what follows we show that if an intrinsic charm component is present in the 
projectile wave function, it will dominate the D-meson production at forward 
rapidities, being more important than the other possible mechanisms. The 
paper is structured as following. In the next section (Sec. |2]) we present a brief 
review of the intrinsic charm models and discuss the predictions of the CTEQ 
collaboration for the x dependence of the charm distribution. In Sec. [3] we 
discuss charm production at forward rapidities and propose a new mechanism 
based on the CGC formalism. In Sec. H] our predictions for the transverse 
momentum distribution of D-mesons produced in pp and p{d)A collisions at 
RHIC and LHC energies are presented. Finally, in Sec. |5] we summarize our 
main conclusions. 
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2 Intrinsic Charm Models 



The charm content of the nucleon sea comes mostly from the DGLAP evo- 
lution of the initial gluon distribution. This process is well understood in 
perturbative QCD. However, as pointed out long ago in Ref. P] (see also Ref. 
[9]) there may be another, higher twist, component which is typical of the 
hadron where it is created. This component is called intrinsic charm (IC). 
A comprehensive review of the main characteristics of the IC models can be 
found in [llj. A well known model was proposed in Ref. [10]. In this model, the 
creation of the cc pair was studied in detail. It was assumed that the nucleon 
light cone wave function has higher Fock states, the first one being \qqqcc >. 
The probability of finding the nucleon in this configuration was given by the 
inverse of the squared invariant mass of the system. Because of the heavy 
charm mass, this probability as a function of the quark fractional momentum, 
P{x), is very hard, as compared to the one obtained through the DGLAP 
evolution. We denote this model by BHPS. 

A more dynamical perspective is given by the meson cloud model. In this 
model, the nucleon fluctuates into an intermediate state composed by a charmed 
baryon plus a charmed meson [19] . The charm is always confined in one hadron 
and carries the largest part of its momentum. In the hadronic description we 
can use effective lagrangians to compute the charm splitting functions, which 
turn out to favor harder charm quarks than the DGLAP ones. The main dif- 
ference between the BHPS and meson cloud models is that the latter predicts 
that the charm and anticharm distributions are different, which generates a 
difference in the momentum distributions of and D~ [20] . 
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In spite of its successes, the idea of intrinsic charm was left aside for some 
years. Very recently it was considered again as an ingredient in the global 
fit of DIS data performed by the CTEQ collaboration [5j. In this update the 
CTEQ group determined the shape and normalization of the IC distribution 
in the same way as they do for other parton species. In fact they find several 
IC distributions which are compatible with the world data. Apart from the 
already mentioned BHPS and meson cloud models, the CTEQ group has tested 
another model of intrinsic charm, called sea-like IC. It consists basically in 
assuming that at a very low resolution (before the DGLAP evolution) there is 
already some charm in the nucleon, which has a typical sea quark momentum 
distribution with normalization to be fixed by fitting DIS data. 

The different parameterizations are presented in Fig. [1] For comparison we 
also present the no-IC distribution, where the charm content of the nucleon 
sea comes from the DGLAP evolution (exclusive charm component). We can 
see that the intrinsic charm distribution is a factor ten (sea like) to twenty 
(BHPS) larger than the no-IC distribution (the long dashed line) and the 
peaks of the IC distributions lie in the large x (> 0.1) region. This behavior 
is so peculiar that it gives us hope to observe IC experimentally. 

3 D-meson production at forward rapidities in the CGC formalism 

Heavy quark production at high energies is usually described in the collinear or 
A;T-factorization frameworks of QCD. Using the collinear factorization, charm 
production is described in terms of the basic subprocesses of gluon fusion 
{g+g ~^ c+c) and light quark-antiquark fusion {q+q — )■ c+c). The elementary 
cross section computed to leading order (LO) or next-to-leading order (NLO) 
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is folded with the corresponding parton distributions and fragmentation func- 
tions. This is the basic procedure in most of the calculations performed, for 
instance, with the standard code PYTHIA. However, at small-x the collinear 
factorization should be generalized to resum powers of asln(s/g|^), where qt 
is the transverse momentum of the final state and y/s is the center of mass en- 
ergy. This resummation is done in the /cT-factorization framework, where the 
cross section is expressed in terms of unintegrated gluon distributions which 
are determined by the QCD dynamics at small-x (For a review see e.g. plj). In 
Ref. the authors have estimated the heavy quark production using the fcr- 
factorization scheme assuming the dominance of the gluon-gluon interactions 
and the presence of high density effects and predicted a much harder open 
charm spectrum. More detailed studies, considering for instance the break- 
down of the /c^^-factorization, have been performed in Refs. [23|. In general, 
these studies assume the dominance of gluon interactions and this implies a 
strong decrease of the spectrum at large rapidities, since in this regime the 
cross section is dominated by very large Xp and small values of Xt- From light 
hadron production, we know that in this kinematical range the cross section 
is dominated by the interaction of valence quarks of the projectile. In other 
words, the cross section is dependent on the partonic structure of the projec- 
tile at large-x. Therefore, the basic assumption of the previous studies about 
heavy quark production at high energies and forward rapidities is satisfactory 
if an intrinsic heavy quark component in the proton wave function is not con- 
sidered. On the other hand, if this component is present and strongly modifies 
the behavior of the charm parton distribution at large x, as shown in Fig. [H 
one must estimate its contribution to the D-meson production cross section. 

In order to calculate the intrinsic charm contribution we generalize the DGLAP 
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CGC factorization scheme proposed in Ref. [T7], which describes quite well 
the experimental data on light hadron production at forward rapidities in 
pp/dA collisions [15J (See also [IS]). In this approach the projectile (dilute sys- 
tem) evolves according to the linear DGLAP dynamics and the target (dense 
system) is treated using the CGC formalism. Apart from [17] we refer the 
reader to [21] for the first derivation of the main formula used here and to 
[22] and [23] for other applications of this formalism, including open charm 
production at RHIC (without intrinsic charm component). It is important to 
emphasize that, as shown in [T7], this approach correctly reproduces the re- 
sults of the linear regime of QCD and coUinear factorization in the low density 
regime. 

In what follows, we assume the existence of charm quarks with large x inside 
the proton, which interact with the target and subsequently fragment into D- 
mesons. This implies that single-inclusive charm production in hadron-hadron 
(nucleus) processes will be described in the CGC formalism by 

xA>(x,,%t) Cc/o(^,(3^), (1) 

\ X p / \ Xp J 

where Pt-, y and Xp are the transverse momentum, rapidity and the Feynman-x 
of the produced D meson, respectively. The variable Xp denotes the momentum 
fraction of the charm inside the projectile. Moreover, fc/p{xp, Q"^) is the pro- 
jectile charm distribution function, taken from the CTEQ group, and D{z, Q"^) 
the charm fragmentation function into D mesons taken from Ref. ^18j. The 
factorization scale is chosen to be = p\. These quantities evolve according 
to the DGLAP evolution equations and respect the momentum sum-rule. In 
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Eq. ([T]), Afjr(x,pT) is the fundamental representation of the forward dipole 
amphtude in momentum space determined by CGC physics. As explained be- 
low, it is obtained from M_a{x-, kx), the forward dipole amplitude in the adjoint 
representation, which is given by: 



-- - J d\e''^' ''' 1 - exp (^-^{r^Ql{x2)y^'^'''''^' 



(2) 



where is the transverse size of the dipole and Qs is the saturation scale given 
by Qs{x) = 1 GeV (— ) ^ with xq — 3 x 10"'^ and A ~ 0.3. For nuclear targets 



contains an additional factor A^^'^. The function 7 in ([2]) is the anomalous 
dimension, which encodes the physical content of the dipole amplitude. We 
follow [I5] and parametrize it as: 

(3) 



-f{w) = 71 + (1 -ii)r^ 



l) + b 



where w = qt/Qs{x2) and 71 = 7(10 = 1). The two free parameters a and 
h were fitted to the RHIC data [U]. The corresponding expression Mp for 
quarks is obtained from Ma by the replacement {r^Ql)"' — > {{C p I C A)'rlQl)'^ , 
with Cf/Ca = 4/9. No K factor (to account for higher order corrections) is 
introduced. The basic characteristic of this model for the forward dipole am- 
plitude is that it explicitly satisfies the geometric scaling property. Besides, in 
this model the large pt limit, 7 — )■ 1, is approached much faster in comparison 
with other phenomenological parameterizations, which implies different pre- 
dictions for the large Pt slope of the hadron and photon yield (For a detailed 
discussion see Ref. [IS]). Finally, we would like to emphasize that this satu- 
ration model is able to describe the experimental data for the light hadron 
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production at mid- and forward rapidities. 



A comment is in order. In a full calculation of D production in the CGC for- 
malism we should consider the IC contribution estimated in this paper as well 
as the contribution which comes from gluon-gluon interactions, as estimated 
for instance in Refs. [2^123] . However, as these contributions strongly decrease 
at large rapidities, we believe that our main conclusions are not modified by 
the inclusion of these terms. 

4 Results 

As the charm densities and fragmentation functions come directly from [5] 
and [T8], respectively, and the dipole scattering amplitude has its parameters 
already fixed in [15], the calculation of the Eq. ([1]) is straightforward and 
parameter free. In Fig. [2] (upper panel) we show the pt distribution of D 
mesons predicted for proton-proton collisions at RHIC energy {^/s = 200 
GeV) and fixed rapidity {y = 4). One can see that the presence of intrinsic 
charm modifies the magnitude of the cross section, as well as the shape of the 
Pt distribution. The behavior of the different models can be easily understood 
if we remember that: 




(4) 



From this expression we see that, in the pt range shown in the figure, at RHIC 
energies and y = 4, the typical values of Xp are in the range 0.4 ^ Xp < 0.95. 
In this range, all IC models predict a larger charm component than the no- 
IC model as shown in Fig. [H This also implies larger values for the cross 
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section calculated with BHPS and meson cloud IC models. The cross section 
calculated with the sea like IC model is smaller, since it predicts, already 
at X ~ 0.6, a very small charm distribution, quite similar to the no-IC one. 
On the other hand, the BHPS and meson cloud models predict a very large 
IC distribution in this kinematical range, which implies larger values for the 
cross section. As expected from Fig. [H the BHPS prediction for the spectrum 
is larger than the meson cloud one. Our predictions for D-meson spectra at 
LHC energy {^/s = 14000 GeV) and fixed rapidity {y = 6) are shown in 
Fig. [2] (lower panel). In this case, the typical values of Xp are in the range: 
0.04 ^ Xp ^ 0.6, which is almost complementary to the RHIC range. In this 
range the sea like model dominates at small Xp, while the BHPS and meson 
cloud dominate at large Xp. These features are directly reflected in Fig. [2] 
(lower panel). 

In order to quantify the magnitude of the enhancement associated to the 
presence of intrinsic charm, we have estimated the ratio between the spectra 
predicted by the different IC models and the prediction without the IC com- 
ponent (no IC). In Fig. Owe present our results for this ratio at RHIC (upper 
panel) and LHC (lower panel), respectively. At RHIC, the enhancement pre- 
dicted by the BHPS and the meson cloud models is large in the whole pt 
range. On the other hand, at LHC the enhancement grows with pt and can be 
of a factor 10 at large values of transverse momenta. The main advantage of 
using formula ([1]) is that is able to account for the non-linear physics, which 
is expected to be more visible with nuclear targets. Therefore we shall extend 
our calculations to forward {y = 3.2) charm production at RHIC in d-Au colli- 
sions. This is shown in Fig. |H Comparing with the proton-proton case, we see 
that in d-Au there is a similar enhancement in the pt spectra and it happens 
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at lower rapidities, becoming easier to be observed. 

So far we have compared IC with no IC predictions, with all the distributions 
coming from the same set (CTEQ) used as input for the same formula ([1]). 
This is the correct procedure to establish the order of magnitude of the ef- 
fect that we are investigating. However, there are other D meson production 
mechanisms that contribute to the experimental background, above which we 
are looking for IC. As mentioned above, charm production at high energies 
is described by gluon fusion and light quark-antiquark fusion. The elemen- 
tary cross section is folded with the corresponding parton distributions. This 
is the basic procedure in most of the calculations performed with PYTHIA. 
This is the dominant production mechanism in the central rapidity region. In 
the forward region it underpredicts the experimental distributions. In Fig. [5] 
we compare our spectra obtained with ([1]) and the CTEQ pdf's for the 
BHPS model (solid line) and no-IC (dashed line) with the predictions of (LO) 
PYTHIA 8.108, supplemented with MRSTMCal (LO) pdf's (dot-dot-dashed 
curve) |25]. We have chosen this particular implementation of PYTHIA be- 
cause it was tuned to reproduce charm production in the central region mea- 
sured by STAR. It should be mentioned that the MRSTMCal gluon distribu- 
tion is a factor two larger than the CTEQ one in the same range of x and 
and therefore the PYTHIA curves presented in Fig. O in dashed lines, should 
be regarded as upper limits. We can observe that for y = 3 (upper panel) 
the IC curve is comparable with PYTHIA prediction at lower values of pt- 
For larger values of pt the IC curve is almost one order of magnitude larger 
than the PYTHIA estimate. For y = 4 (lower panel) the IC distribution is al- 
most a factor 10 larger than the PYTHIA one in all the pt range considered. 
Of course, at this stage there are still many uncertainties in this comparison 
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and no definite conclusion can yet be drawn from Fig. [51 However our results 
are encouraging. They strongly suggest that the future measurement of open 
charm forward pt spectra may show an enhancement over standard PYTHIA 
predictions, which is very hard (if at all!) to explain without intrinsic charm. 

Some comments are in order here. Firstly, an alternative probe of the IC may 
be the study of the pr dependence of the ratio (-D/vr) between the D and vr 
production cross sections at a fixed rapidity. We expect an enhancement of 
this ratio, as well as a different pt behavior, if intrinsic charm is present. Sec- 
ondly, although we have only considered positive rapidities, the pr spectra are 
identical for negative rapidities in the case of pp collisions, i.e. our predictions 
depend only of \y\. However, for pA collisions, this symmetry is broken, due 
to the different magnitudes of the high density effects in the proton and in 
the nuclei at the same x. As it is well known, these effects are amplified by 
nuclei [1]. Interestingly, the study of the D meson production in the nucleus 
fragmentation region can be useful to determine the presence of IC and its 
behavior in a nuclei. Finally, it is important to emphasize that the study of 
D-meson production at forward rapidities as a probe of the intrinsic charm 
component was proposed and studied many years ago in Refs. piGpO] . How- 
ever, these authors considered that the intrinsic charm is put on shell and 
emerge in a real charmed particle through soft non-perturbative interactions. 
The charmed meson in the final state is generated by the coalescence of the 
intrinsic charm with co-moving light spectator quarks. Basically, the process 
of D-meson production at forward rapidities, associated to the presence of 
intrinsic charm component, proposed in these references is essentially non- 
perturbative. In contrast, the mechanism proposed for the first time in this 
paper is based on a perturbative treatment of the interaction of the intrinsic 
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charm in the projectile with the Color Glass Condensate in the target. 

During the revision of this work, a closely related paper appeared [26j , propos- 
ing the search for IC effects in distributions of D mesons produced in proton 
- antiproton collisions at Tevatron and measured by the CDF collaboration. 
In this kinematical set-up one can also have access to the large x component 
of the projectile wave function. Indeed, recalling we observe that large Xp 
is dominant in the production of D's either with large y and moderate pt, 
which is the region explored here, or with large px and low y, which is the 
region considered in Both works reach similar conclusions and may be 
regarded as complementary. 



5 Summary 



Although the direct measurements of heavy flavors in DIS are consistent with 
a perturbative origin, these experiments are not sensitive to heavy quarks at 
large x. Therefore, it is fundamental to study other observables which may be 
used to determine the presence (or not) of an intrinsic heavy quark component 
in the hadron wave function. In this paper we have proposed a new mechanism 
for the production of D mesons at forward rapidities which assumes that the 
appropriate formalism in this kinematical range is based on CGC physics, 
and demonstrate that the resulting transverse momentum spectra are strongly 
dependent on the behavior of the charm distribution at large Bjorken x. Our 
results indicate that the intrinsic charm hypothesis can be tested in pp {p{d)A) 
collisions at RHIC and LHC. 
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Fig. 1. Intrinsic charm momentum distribution in the proton with different IC mod- 
els taken from [5]. 
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Fig. 2. Transverse momentum distribution of D mesons produced in p-p collisions 
at RHIC and LHC. The spectra are calculated with ([T]) and the four charm densities 
taken from [3] and shown in Fig. [TJ 
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Fig. 3. Ratio between the different spectra calculated assuming IC and the no-IC 
spectrum shown in Fig. [2l 
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Fig. 4. Upper panel: spectra calculated for deuteron-gold collisions with several IC 
models and with no IC. Lower panel: ratio IC/ no IC. 
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Fig. 5. Comparison between our spectra calculated assuming IC and the predictions 
of PYTHIA. 
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